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AN INVESTIGATION OF THE EFFECT OF TWIST ON THE COUPLED 



BENDING-TORSION VIBRATIONS OF A CANTILEVER BEAM 



SUMMARY 

An experimental study was made of the effect of twist on the natural 
frequencies, coupled in bending and torsion, of a cantilever beam. The shape of 
the beam was selected to provide a high degree of coupling. An analytic procedure 
for determining these frequencies was also developed and is included in this report. 

It was determined that coupling without twist reduced all frequencies from 
their uncoupled values. The fundamental frequency was found to remain relatively 
constant as the beam was twisted. All natural frequencies higher than the 
fundamental were lowered as the total twist was increased to 15 degrees, but 
remained relatively constant with further twist. 



INTRODUCTION 



With the advance of helicopters and turbo-machinery, the problem of twist 
and its effect upon the natural vibrations of a beam becomes of more and more 
interest. Much work has been done to determine the effect of twist on beams 
having no coupling between torsion and bending, and the effect of a high degree of 
bending-torsion coupling on untwisted beams. This report shows primarily the 
effect of twist on the natural frequencies of a particular beam highly coupled in 
bending and torsion. It also shows the node lines for all the natural frequencies, 
and indicates the effect of bending -torsion coupling on the natural frequencies in 
the untwisted beam. 

A beam with dimensions as shown in Fig. 1 is considered. This particular 
shape was chosen to give a maximum amount of bending -torsion coupling, 
relatively easy tooling problems, and readily obtainable sectional and elastic 
properties. In this report, the term coupling refers to the interaction of the 
torsional, flapwise bending, and chordwise bending types of vibrations. It was 
planned to use both experimental tests and analytic means to determine the 
natural frequencies resulting from the combination of twist and coupling. Only 
experimental results however, are contained in this report. Difficulties en- 
countered in preparing the digital computer program prevented the inclusion of 
analytic results. The mathematical method used is outlined in the Appendix, 

This method is based on one being used in a report of Isakson and Eisley 
currently being prepared for the University of Michigan Research Institute. It 
involves the introduction of coupling terms from the equations of motion developed 
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by Houbolt and Brooks, Ref. 1, into the matrix method of solution of Targoff, Ref. 2 &3. 
The Targoff method was described and used in a report of Isakson and Eisley, Ref. 4. 

The coordinate system used is shown in Fig* 2. The flexural center of the 
cross section was located according to a formula from Roark. Ref. 5, which formula 
was developed from Ref. 7. The elastic axis of the beam was maintained as a 
straight line by mechanically twisting the beam about this axis. The nomenclature 
and sign convention for the displacements, shears and bending moments used in the 
analysis are shown in Fig. 2. All parameters used or needed are listed in the table 
of symbols. The non-dimensional forms of these parameters used in the analysis 
are also listed in the table of symbols. The mass and stiffness distributions are 
uniform for the length of the beam. 

In gathering the experimental data, tests were made at total angles of twist 
of 0°, 5.5°, 7.90, 12.9°, 20.7°, 2G.7°, 30.4°, and 40. 3°. 

In conducting this investigation, the assistance and guidance of Professors 
Eisley and Isakson, of the Department of Aeronautical Engineering of the University 
of Michigan, were of great value. 
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SYMBOLS 



A 

A 

B i 

B„ 



GJ + EB X X ( j3' ) 2 

A_ 

El, 



cross section constant 



/’” + t -l)Ji 

cross section constant f f] (Y\ + ~ k* ) d t} 



EB R’ 
2 



E V EI 2 



G 

GJ 



V 1 * 



In . I 



1 ’ l s 



l S 



eb 2 b ;/ 

ei x /a 

Young’s Modulus 

bending stiffness about major and minor 
principal axis of cross section respectively 
Shear Modulus 

St Venant’s torsional stiffness 
moments of inertia about major and minor neutral 
axes, respectively (both pass through centroid of 
cross sectional area effective in carrying tensions). 
Moment of inertia of cross section about elastic axis 
mass moment of inertia per unit length about 
and J axes respectively 

Li 

/°. R 



h . 

AR' 

Torsional Stiffness constant J r 



Vi, 



1 + 



n J, 



(Ac? *) (CA o < A ) 



(Ref t) 
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Bending moment in flapwise direction 
Bending moment in chordwise direction 
Torque about elastic axis at any cross seotion 
Beam length 

Shear load in flapwise and chordwise direction 
respectively 



mode constant for determination of frequency 
distance between area centroid of tensile member 
and elastic axis, positive for centroid forward 




acceleration due to gravity 

polar radius of gyration of oross-sectional area, 
effective in carrying out tensile stresses, about 
elastic axis 




El, 

LiM l 

Uit 



k 



A 



R 



length of blade segment 



1 _ 

R 



mass of blade segment 



^ j 


total loadings per unit length in x l ij , and ^ 
directions 


Qx 


total torque loadings per unit length in x , f\ , (Mid £ 
directions 


t 


thickness of cross section at any chordwise station 


X 


coordinate in direction of X axis measured from 
the root along the elastic axis 


( 3 


blade angle between major axis of cross-section 
and the "y" axis for <fi - 0 , 
positive when leading edge is up. 


(J 


rate of twist of blade per unit length 


A/3 

Y 


increment in Q between blade segments 

■Jei: 

Uli 


°y > <b 


displacement of elastic axis in "y" and " z " direction, 
respectively 


0 / i Pi. 


displacement of elastic axis in direction of minor and 


e 


major principal axes of the cross-section, respectively 
total twist in blade between X - 0 and X = R 


f-i 


cross-sectional coordinate; Y] axis lies along major axis; 
£ axis is perpendicular to major axis and passes through 


X 1 

P 


elastic axis. 

UlV° ft 
El. 

angle of twisting deformation, positive when leading edge is up 


p 


mass per unit length of blade 


o 


z 



' A 

Note: all zero subscripts refer to root section 
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Parameter Values for Beam Under Consideration 



Values Independent of Twist 







6 


B 1 


93.941 


in 


B 2 


-11.50 


. 5 

in 


E 


10. 5 x 10 6 


lb/in 2 


El i 


.7959 x 10 6 


lb in 2 


EI 2 


32. 31 x 10 6 


lb in 2 


G 


4.0 x 10 6 


lb/in 2 


GJ 


/•» 

1.1864 x 10 


lb in 2 




. 0758 


. 4 

in 


<2 


3.125 


. 4 
in 




15.552 


. 4 

in 




. 18723 x 10~ 4 


2 , 2 
lb sec /in 




38.4099 x 10' 4 


IK 2 /- 2 

lb sec /in 


J 


.2966 


. 4 
in 


R 


26.0 


in 


^A 


2.350 


in 


g 


32.2 


ft/sec 2 


k A 


2.635 


in 


1 


2.6 


in 


m 


.00145 


2 

lb sec /in 


P 


. 00055575 


2 . 2 
lb sec /in 
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Values which are dependent on twist 



Total twist: 


0° 


15° 


o 

30 


A (lb in^) 


1.1864 x 10 


1.2699 x 10 6 


1. 5196 x 10 


C (radians) 


0 


874780 


-1.396937 


6 (degrees) 


0 


15 


30 


6 (radians) 


0 


.2388 


.4776 


[\j3 (radians) 


0 


. 02388 


. 04776 


(3 (rad/in) 


0 


. 00920 


. 01838 



7 



Non Dimensional Parameters 






if 



^A 



1 

? 

For 9 - 



A 

c 

X 

Y 



Values independent of twist 

.498356 xlO* 

102.2392 x 10 4 

. 09038 
.1013 
.1 

.0243 

1 



Values dependent on twist 



1.4906 



o 



0 

0 



1.59550 
-. 374780 
. 12617 
. 0472554 



1.90440 
-1.39694 
. 50469 
. 192792 
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EQUIPMENT AND PROCEDURE 



The model used in the present investigation was made of 2024 T 4 Aluminum, 
machined from 6" x 1" stock. It was machined to within . 015 inch of the given 
beam dimensions. 

Four inches at the end of the model was left rectangular for mounting purposes 
(Fig. 1). 

The bar was mounted with four bolts between two modified channel sections, 
which in turn were fastened to a double box section bolted to a structural member 
in the wall of the building. The double box section was made by welding 1/2” inch 
steel plates between the outer edges of the flanges of a 6 x 6, 5/8 inch steel I beam. 
(Figs. 3 and 4). The channel sections were later reinforced by welding a 0. 5 inch 
steel gusset as shown in Fig. 3. This was done to observe the effect upon the 
natural frequencies of substantially stiffening the mounting. No measurable change 
in frequencies was noted. This provided some indication that the mounting was 
sufficiently rigid to give results that closely approached a true cantilever. 

MB Vibration Test Equipment, Model T1 - 32034 was used to vibrate the 
beam. This machine has a frequency range of from two to 20, 000 cycles per second. 
The vibrator itself was mounted to the floor and fastened to the beam with an 
aluminum tube. This tube was fastened to the vibrator with a flexure joint and to 
the beam with a universal joint (Figs. 3 and 4). 

Two universal joints were used originally, one at each end of the aluminum 
tube, but this arrangement allowed too much slack in the connection. Two flexure 
joints were also tried in the same manner as the universals, but this arrangement 
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was too stiff to permit the beam sufficient freedom in torsion. 

To obtain the desired twist, special fittings were made to fasten to each end 
of the bar so that it could be twisted about its elastic axis. A Riehle torsion testing 
maohine of 120, 000 in. lbs. capacity was used to twist the beam. The twisting rig 
as shown in Figs. 5 and 6 was used to transmit the torque. 

The method of determining natural frequencies consisted of the following steps: 

1) A frequency was approximately located by listening for and observing the 
amplitude maximum. In the case of the fundamental frequency this method provided 

an accurate answer by noting the frequency at which the measured maximum amplitude 
of the free end of the beam occurred. 

2) The beam had been designed so that uncoupled natural frequencies of torsion 
and chordwise bending occurred in the vicinity of a natural uncoupled flapwise 
bending frequency. When the beam was being vibrated in the vicinity of a band of 
natural frequencies, fine sand of silicone flour was scattered over the surface of the 
beam. The frequency was then adjusted until sharp node lines appeared and concided 
with maximum amplitudes of vibration. The maximum amplitudes were recognized 
by first turning down the amplitude control of the vibrator until the sand or powder 
was barely agitated, then adjusting the frequency until the maximum agitation of the 
particles appeared. This method permitted the determination of amplitude maximal - 
in a band of frequencies where sharp node lines did not disappear, but only shifted 
their positions. Where several maxima u appeared in a band, the strongest was 
assumed to be the frequency of the flapwise bending mode. This assumption was 
made considering that a flapwise node was designed to occur in eaoh band, and that 
due to the shape of the cross section, the flapwise vibration would be the dominant one. 
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It is the frequency change due to coupling effects as influenced by twist 
rather than the absolute value of the frequencies that is important in this report 
Several factors occur in the experimental procedure which would tend to produce 
an experimental result different from one arrived at by calculations. Throughout 
the experimental runs, attention was directed toward keeping these factors constant. 

The connecting rig between the beam and the shaker added end mass effect; 
this is discussed in the next section of this report. It was noted that varying the 
tightness of the bolts at the mounting even slightly could produce a slight variation 
in frequency. Of course, for all runs the beam was mounted as securely as possible 
so that the net effect of variations from one test to another was very small. 

For the evaluation fo the analytic part of this report the Royal Precision 
Electronic Computer LGP-30, manufactured by the Royal McBee Corporation was 
to be used. A matrix program was set up as shown in the Appendix. The computer 
evidently oould not handle the problem as it was programmed. A great deal of 
time had been spent on this part of the investigation and Insufficient time remained 
to pursue it further. 
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RESULTS AND DISCUSSION 



The effects of twist upon the natural frequencies of the beam are shown in 
Figs. 7 through 12. In Figs. 7 through 9, the absolute values of the natural 
frequencies are plotted versus the angle of twist. In Figs. 10 through 12, instead 
of absolute frequencies, the ratio of each frequency at any degree of twist to the 
corresponding value for the untwisted beam is plotted versus the angle of twist. 

The frequency ratios are grouped according to the assumed predominant motion, 
i.e. flapwise, chordwise or torsional. 

Table I lists the uncoupled natural frequencies of the beam calculated on the 
basis of geometry and mass of the beam. These are compared with the corresponding 
natural frequencies observed when vibrating the untwisted beam. A coupled rrasonant 
frequency appears in the frequency spectrum corresponding to each calculated un- 
coupled natural frequency. It was assumed that the type of motion associated with 
each uncoupled natural frequency, remained a6 the predominant motion in the 
corresponding coupled natural frequency of the untwisted beam, and would continue 
to predominate as each frequency changed with twist. 

The discussion of the effect of twist is based primarily on Figs. 10 through 12, 
the frequency ratio plots. In this discussion each frequency is identified by its 
predominant motion. 

Twist has very little effect upon the fundamental frequency (first flapwise 
bending) of the beam. This frequency increases very slightly as the twist is increased 
through 30 degrees. The data point obtained at 40.4 degrees of twist indicates that 
there is a possibility of a more rapid rise in frequency as the twist is increased 
beyond this point. 
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The second flapwise frequency drops steadily as twist is increased to 
approximately 15 degrees; it remains constant as twist is increased further. 

The third and fourth flapwise frequencies decrease as twist is increased 
to 10 degrees. For the next 10 degrees of twist, these frequencies rise slightly, 
then remain constant as twist is increased further. 

The first torsional frequency, corresponding to the fourth natural coupled 
frequency of the beam, drops off rather rapidly as twist is increased to 15 degrees, 
it then decreases slowly as twist is further increased. The second torsional 
frequency decreases until about 15 degrees of twist is obtained, then tends to 
remain constant. 

The third torsional frequency decreased to 15 degrees of twist then appeared 
to rise as twist was further increased. Data for the third torsional frequency was 
considered less reliable than the other data, because of the difficulty encountered 
in obtaining its resonance point and in obtaining sharp node lines. 

Only the first chordwise frequency, corresponding to the third natural 
coupled frequency of the beam, was obtained. It dropped sharply until a twist of 
about 20 degrees was obtained; it then appeared to remain constant or rise slightly 
as further twist was introduced. 

The shape and location of the node lines for the different frequencies are 
shown in Figs. 13 through 20. These were obtained when determining the natural 
frequencies. Although in some cases a band of frequencies occurred in which the 
node lines remained sharply defined and only changed their orientation and position 
with frequency, resonant peaks within the band were strong enough to permit 
determination of resonant frequencies. The occurence of frequency bands is 
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indicated in Figs. 13 through 20. 



The shape and location of the node lines was significant in that they helped 
in the identification of the type of motion. It was noted that the node lines tend to 
move toward the free end of the beam as twist is increased to about 20 degrees. 

With further twist they remain relatively constant or even move slightly away 
from the free end. The position of the node lines in each group of frequencies 
reflects the flapwise bending mode occurring in that group. 

The node line illustrations show the appearance of what was assumed to be 
resonant frequencies of plate bending. This assumption was made since no 
natural frequencies based on flapwise bending, chordwise bending, or torsion, were 
predicted to oocur in this range. This plate bending first appeared near 300 cycles 
per second and became more pronounced as twist was increased. At high angles of 
twist, (20 to 30 degrees), it began at approximately 276 cycles and produced 
distinct node lines through a band of 400 oycles per second. Twist had little effect 
on the mean frequency of this mode as can be seen in Fig. 21. 

This plate bending appeared again at about 680 cycles for a twist of 25.7 degrees, 
becoming evry distinct at a twist of 30.4 degrees. Another plate bending mode was 
seen to oocur at 960 cycles and first appeared at a twist of 5. 5 degrees. 

Fortunately, the plate bending resonant frequencies appeared between the 
other bands of frequencies and though they undoubtedly influenced the shape of the 
node lines in the upper range of frequencies, it is felt that they did not greatly 
hinder the accurate determination of resonant frequencies in modes of vibration 
which are of interest in this paper. 

In shaking the beam, oertain frequencies above the fourth natural frequency 
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did not give distinct nodes. In these cases, the natural frequencies were plotted 
as a band rather than as a distintc point (Figs. 7 through 9). This difficulty was 
due not only to the proximity of the flapwise, chordwise, and torsional uncoupled 
frequencies, but also to the appearance of plate bending which was discussed above. 
Other inacouracies were due to some unavoidable looseness in the vibrating 
mechanism, primarily the universal joint. This looseness also limited the frequency 
range obtainable, since the vibrating movement was lost in transmittal at high 
frequencies. The tube and fittings which connected the bean to the vibrator also 
added a slight end mass effect to the beam. By adding additional known mass 
weights to the end of the beam, it was determined that the effect of the vibrating rig 
reduced the fundamental frequency by 1.6 cycles per second. The data taken was 
not corrected for thi3 effect, however, since the effect at higher frequencies was 
not known. The frequencies might also have been slightly low due to the impossibility 
of obtaining a perfectly rigid support at the root. The base of the beam however, 
was considered sufficiently rigid so that this effect was minor. 
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CONCLUSIONS 



1. In all cases, the effect of bending-torsion coupling, independent of twist, 
appears to reduce the natural frequencies of the beam. 

2. The fundamental frequency is relatively insensitive to twist, demonstrating 
only a slight rise as the angle of twist is increased. 

3. All natural frequencies above the fundamental one are lowered as twist is 
increased to approximately 15 degrees. 

4. Beyond 15 degrees of twist the natural frequencies are relatively independent 
of further twist. 



I 
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APPENDIX 



METHOD OF ANALYSIS 

The analysis is basically the one which is outlined in Ref. 3 and the 
Appendix of Ref. 4. It is extended to include torsional vibration and bending- 
torsion coupling and is applied to the non-rotating case. 

The beam is divided spanwise into ten equal segments. The mass of each 
segment is assumed concentrated at its center, and the bending stiffnesses EI^ 
and EI 2 and angle of incidence yS are assumed constant between masses. The 
twist of the beam is accounted for by relative rotations of adjacent uniform bays 
(between masses) about a spanwise axis, the change in angle, A ft , being equal to 
the total twist in a segment and occuring just outboard of the mass. 

With torsional motion included, the following quantities , at any point along 
the beam, may be expressed as a column matrix: 




S* 

*2 



Q 

< 0 > 



1 



The elements of this matrix vary from one station to the next along the beam 
in the following fashion: | A j = [ R ][ F ][ F ]( A }n where [r] [e] [f ] 

are 10 x 10 matrices representing linear relationships, described below, between the 
elements of anc * corresponding elements in |aJ ^ ( . Coupling effects are 

introduced into these relationships by consideration of the differential equations of 
coupled motion developed by Houbolt and Brooks in Ref. 1. 

In the jj^J matrix whioh relates the |a| matrices across a weightless 
section, the elements are obtained as follows: 



(a) 


w (ni-i) _ 

v i 


v, w 




(b) 


jvj Com) _ 


M, (n) + V, w t 






From Ref. 


1, at any station: 








E I, (-v"sm/3 + w 


"cos/3) 






= El .(zsip + s; 


- 






0 

II 

<=CL. 

»— T 

LU 

1 1 


for eac h 




( E I, 


h) M (n) w (n) 

= n, + v, * 




(c) 


w(n-n) 

O/ 


s; ln) - 


Art - 
tr; 


(d) 


= 


s, M - s; ln) i + 


nA 1 . 

OCT 


(e) 


\4 (n,,) = 




fc" U A | 


(f) 


ri i !nH) = 


ri 2 (n) + v 1 (n, t 






From Ref. 


1, at any station: 






= 


E (v (, cos|3 r w" in 


n(3) -EB. 



v/ 0 V 

zk i, 
v, ln> t J 






= E I, E BjP 

= E IS" - tBiP'0' 



2 



For the untwisted section ^3 » 0 , but the second term in the above equation 
is not dropped since it represents the ooupling between chordwise bending 
and twist which can only appear at this piont in the relationships. For this 
term, j3 must equal the amount of twist that actually exists in a segment 
having length x. 

(E 1^“ - : Mr’ 1 ’ V 2 <n)x a) 

<g) Q (nM > = Q (n) 

From Ref. 1, at any station: 

Q =L&J+ e <P' - E + w'iin/3) 

= [ G J *■ E 6, ( p') z ] (f)‘ - E5i S" (p-0 U, each S tat/on) 

The {I appears in a coupling term and is retained as before. 

[gj 0'- eb^'s; = Q (n> ^ 

Equations (1) and (2) are solved simultaneously for and 0 ; 

from & L ' , Si and are obtained by integration. 



( 2 ) 



<b) 



/ (dm) 



(i) 






(Dfl) 



(A I1 ln, i- EBj /3'Q (n) ) £ + A vf’f- 

EU [A- i£M:> 1 ] 

Eh 



+ 



Si 



M 



K (n > + *1 



ErJ a- lHi£l l j 

* EI, 1 



(n) 



(J) 



<b 



(n*i) _ 



0 (n) f 



t- EBj ft' / 

Eh l M < 



W o ^ Vi n) 2 z 



e t 



) 



A - iiMl i 

Eh 



0 



W 






3 



where A = GJ + EB^(|# )* 

Equations (a) through (j) provide the elements of the £ eJ matrix; these 
elements are put into non-dimensional form and the result is the fOxlO 
[ e] matrix whose terms are all zero except for the following: 




The quantities B and B are those defined by Houbolt and Brooks, Ref. 1, 

JL u 



and were evaluated for the particular beam. 

The elements of the matrix which relate the | a| matrices on either 



side of a concentrated mass are determined as follows: 
(n»<) 



V,' 



(a) 



= v, 1 "’ 




m ^ m Qft 0 


= v, (n) 


+ 


mu> 2 S"/ n * + m id 






m k U) 2 " £ Si 



,(n) 



(b) 



The second term is one resulting from the entire segment's undergoing a 
change of slope when the beam bends. 

k^ = radius of gyration of a segment about a chordwise axis; m = mass of a 



segment. 






(c) 


C / (0f») 


c ' (n) 




b t 


o, 




<d) 




6, ln) 






V i ‘ n '" ) = 


v 2 tn) + 


r (n) 
m 


(e) 




\/i n) H- 






a lnM) = 


M/" 1 + 


m kx 


(f) 


— 


Mi" ) + 


m 1 $ 2 . 



The second term is a rotary inertia term as before; 



k - radius of gyration of a segment about a vertical axis; m=.mass of a 

Z 



z 






segment. 






r-i(n-l) 




r* ' ( n ) 


(g) 0 Z 




*1 


« S * 1 "' 0 


— 


S L (n) 


(i) Q (n "> 




Q ^ + m -+• m kj 0^ 




— 








where k is the polar radius of gyration of 

u 






the segment. 


0) 0 £ ”'° 


— 




Equations (a) through (j) are put into non-dimensional form and provide the 


elements of the 


M 


matrix: 
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1 



F = F = F -F — F =F — F =F =F = F 
11 22 33 44 55 66 77 88 99 1010 



F =• F = f P V 
14 58 r 1 



Z So > i 



F =F = 3 i. <=« 

110 94 r * 

p _ p — 4 f 

F 23 F 67 >12 



F s>io - (I T fI 5 ) f ** 

The [R] matrix serves to rotate the coordinate axes through the angle Zi /3 
and appears as follows: 



R 

R 

R 

R 



11 

15 

51 

99 



“ R 22* R 33 
= R 26 = R 37 
” R 62= R 73 



=• R 



1010 



= R 
= R 
— R 
1 



44 

48 

84 



= R * R = R 
55 66 77 

= - sin a/3 

= sin /\f3 





The analysis then continues in the same manner as in Ref. 4, exoept that the 
matrices are now 10x10 and yield a 10x5 product when the boundary conditions at 
the tip are introduced. 

The final results are obtained through evaluation of a 5x5 determinate. 
Introduction of trial values of 7^ ^ L 7* “ f ( u> ) ] , produce points 

which define the residual curve. Zero values of this curve define the natural 
frequencies. 
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TABLE I 



COMPARISON OF COUPLED AND UNCOUPLED FREQUENCIES 



F requency 
(Cycles per second) 


Type of Vibration 


Uncoupled* 


Coupled 


30.6 


25.9 


1st Flapwise 


169 


155 


1st Torsion 


195 


163 


2nd Flapwise 


198 


177 


1st Chordwise 


507 


488 


2nd Torsion 


536 


518 


3rd Flapwise 


844 


— 


3rd Torsion 


1054 


925 


4th Flapwise 



•Does not include the effect of end masses. 
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Effect of Twist on Natural Frequencies 
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Fig. 13 
Node Lines 
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Fund.Freq. =• 25.9 
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Fig. 14 
Node Lines 
6 - 5 . 5 ° 



Fund. Freq.=26.0 
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Fig. 15 
Node Lines 
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9 = 7.9° 

Fund. Freq. -26.1 



Node at end of beam 
weakened at 460 cycles 



Node at center of beam 
weakened at 480 cycles 
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Fig. 16 
Node Line a 
0 =12.9° 

Fund. Freq. = 25.8 
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FIs. 17 
Node Lines 
9 = 20.7° 

Fund. Freq.=26.0 
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Fig. 18 
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Node Line a 

0 « 26.7° 

Fund. Freq.=25.7 
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Fig. 19 
Node Lines 
6 =30.4° 

Fund. Freq,=26. 8 
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Fig. 20 
Node Lines 
0 = 40.8° 

Fund. Freq. =27.1 
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